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ABSTRACT 

We explore the stability of different galaxy light concentration indices as a function of the 
outermost observed galaxy radius. With a series of analytical light-profile models, we show math- 
ematically how varying the radial extent to which one measures a galaxy's light can strongly affect 
the derived galaxy concentration. The "mean concentration index" , often used for parametrizing 
high-redshift galaxies, is shown to be horribly unstable, even when modelling one-component 
systems such as Elliptical, dwarf Elliptical, and pure exponential disk galaxies. The C31 concen- 
tration index performs considerably better but is also heavily dependent on the radial extent, and 
hence exposure depth, of any given galaxy. We show that the recently defined central concentra- 
tion index is remarkably stable against changes to the outer radius and observational errors, and 
provides both a meaningful and reliable estimate of galaxy concentration. The Sersic index n 
from the r^/" models is shown to be monotonically related with the central concentration of light, 
giving the index n a second and perhaps more tangible meaning. With a sample of Elliptical and 
dwarf Elliptical galaxies we present correlations between the central light concentration and the 
global parameters: luminosity (Pearson's r — —0.82), effective radius (r = 0.67), central surface 
brightness (r = —0.88), and velocity dispersion (r — 0.80). The more massive Elliptical galaxies 
are shown to be more centrally concentrated. We speculate that the physical mechanism behind 
the recently observed correlation between the central velocity dispersion (mass) of a galaxy and 
the mass of it's central supermassive blackhole may be connected with the central galaxy con- 
centration. That is, we hypothesize that it may not simply be the amount of mass in a galaxy, 
but rather how that mass is distributed, which controls the mass of the central blackhole. 

Subject headings: black hole physics, galaxies: fundamental parameters, galaxies: kinematics and dy- 
namics, galaxies: nuclei, galaxies: photometry, galaxies: structure 

troduction of the more commonly used C31 index 
(see, for e.g., Kent 1985; Gavazzi, Garilli, & Boselli 
1990; Moriondo et al. 1999). 

Okamura, Kodaira, & Watanabe (1984) ex- 
plored other fractional ratios and also introduced 
what they called the "mean concentration index" . 
This was a ratio of the luminous flux enclosed by 
two different radii rather than a ratio of radii (see 
section 2.2). Using an outer isophotal level of 26 
V-ma.g arcsec"^, Okamura et al. described how 
these indices are dependent on the intrinsic (or 



1. Introduction 

Estimates to the central concentration of galaxy 
light for the parametrization of galaxies go back 
many decades (Morgan 1958, 1959, 1962; sub- 
sequently known as the Yerkes system). Eraser 
(1972) made this method of classification quanti- 
tative with the introduction of the concentration 
indices C21 and C32, defined as the ratio of radii 
that contain 50% and 25%, and 75% and 50%, 
of the total (asymptotic) galaxy luminosity, de 
Vaucouleurs (1977) expanded on this with the in- 
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mean) surface brightness of the individual galax- 
ies. At least for Elliptical galaxies, or rather, one- 
component galaxies or bulges - which is what we 
wish to explore hero this can be directly trans- 
lated into a dependency on the number of effec- 
tive radii one samples. It is desirable to perform 
this translation because one replaces a dependency 
on two quantities, namely, the intrinsic surface 
brightness of each galaxy and the faintest surface 
brightness level observed, with one quantity (i.e. 
the number of scale-lengths or effective radii ob- 
served). It is this dependency which we wish to 
explore. Ideally, one would like to have a concen- 
tration index which is as independent as possible 
on this quantity (which is dependent on, amongst 
other things, the observational exposure details). 

We use the Sersic (1968) r^" law as a model- 
dependent way to illustrate the various concentra- 
tion indices presented in Section 2. (The concen- 
tration indices, however, can be measured inde- 
pendently of this, and any, model.) All the con- 
centration indices are dependent in some way on 
the extent to which the galaxy radial profile is 
sampled, and this dependency is revealed in Sec- 
tion 3. Truncating analytical light-profile models, 
we perform a comparative study of three different 
concentration indices. Additionally, using a range 
of observed Elliptical galax;y profiles obtained from 
deep exposures, we again explore the stability of 
the concentration indices. 

In Section 4 we investigate the relationship be- 
tween the galaxy light concentration and the other 
physical galaxy properties such as: luminosity, ef- 
fective radius, central surface brightness, and ve- 
locity dispersion. Under the assumption that El- 
liptical galaxies arc homologous systems, that is, 
assuming they all obey the r^/^ law, the concen- 
tration index of every Elliptical galaxy should be 
the same. However, research over the last decade 
has shown that Elliptical galaxies are not homol- 
ogous and we present, in some cases, the first ever 
correlations between galaxy concentration and the 
global properties of Elliptical and dwarf Elliptical 
galaxies. 

2. Concentration indices 

We use the Sersic r^/" law as a mathemati- 
cal means to illustrate the various concentration 
indices. The r^/" luminosity-profile model has 



been shown to provide a good representation to 
the distribution of light in both Elliptical galaxies 
(including the dwarf Ellipticals) and the bulges 
of Spiral galaxies (Caon, Capaccioli, & D'Onofrio 
1993; D'Onofrio, Capaccioh, & Caon 1994; Young 
& Currie 1994; Andredakis, Peletier, & Balcells 
1995), and is also appropriate for describing a Spi- 
ral galaxy's disk. Although it's appearance in the 
literature is becoming more frequent, we summa- 
rize it's form below and introduce the parameters 
we will subsequently use throughout the paper. 

The de Vaucouleurs (1948, 1959) r^/^ radial 
intensity profile /(r) was generalized by Sersic 
(1968) to give the r^/" law where 

/(r) =/(0)exp-*''(''/''-)'^" . (1) 

7(0) is the central intensity and re the effective ra- 
dius enclosing half of the flux from the model light 
profile. The quantity 6„ is a function of the shape 
parameter n - which defines the global curvature 
in the luminosity profile - and is obtained from 
solving the expression T{2n)=2^{2n,bn), where 
r(o) and 7(a, x) are respectively the gamma func- 
tion and the incomplete gamma fimction. The to- 
tal luminosity, Lt, associated with an r-^^" profile 
that extends to infinity can be written as 

^T = /(0)r,%„,|jr(2n), (2) 

where r,. mofj is the effective half-light radius of the 
model. For elliptical (that is, non-spherical) galax- 
ies, ^-irnod = ''e,maj(l-e), whcrc € is the cllipticity 
of the galaxy and re,maj is the major-axis half-light 
radius. The outer radius to which one (reliably) 
measures a galaxy is of course a function of expo- 
sure time, telescope aperture, etc. Denoting this 
outer finite (or final) radius as rfin, the luminosity 
enclosed by this radius is 

"n \ V ' e,mod / / 

(3) 

Having introduced a model which can be used 
to represent the observed range of structural pro- 
file shapes in bulges, we can now proceed to de- 
scribe the degree of concentration of light in these 
systems. The following definitions for the various 
concentration indices can all be applied without 
reference to the above (or any) light-profile model; 



2 



but it is of course insightful to use a parametrized 
model. We give below definitions or procedures 
to obtain the various concentration indices as a 
function of the outermost observed radius. 

2.1. The concentration index C31 

Due to the popular use of C31 over C32 and C21 , 
we will focus on this particular concentration in- 
dex related to the ratio of radii. The C31 index 
can either be obtained directly from the image, 
or using equation 2 and 3 when the observed lu- 
minosity profile is well fitted with an r^/" model. 
Because light-profile models can first be convolved 
with the relevant point spread function before fit- 
ting to an observed profile, the latter approach 
has the strong advantage that corrections for see- 
ing are already taken into account. Dividing equa- 
tion 3 by equation 2 gives the fractional luminos- 
ity, denoted here by x, contained within the out- 
ermost observed radius rfin. One then solves the 
ratio of equation 3 and equation 2 for the (new) 
values of rfin/re,mod that give a fractional luminos- 
ity of 0.75.T and 0.25.x. The results of doing this 
are shown in Section 3. 

2.2. The mean concentration index 

The fundamental parameter in the classifica- 
tion system of Abraham et al. (1994) is basically 
the concentration index defined by Okamura et al. 
(1984) and Doi et al. (1993). ft is a luminos- 
ity ratio between the flux enclosed by some in- 
ner radii and the outer-most radii, and has been 
parametrized in Abraham et al. (1994) such that 



outer radius rfin such that 



C{a) 



(4) 



Here, represents the intensity in the pixel (i, j), 
and E{a) denotes some inner radius which is a 
(0 < a < 1) times the outer radius which has 
been normalized to 1. Following Doi et al. (1993) 
Abraham et al. (1994) used a value of 0.3 for a. 
In what follows in Section 3, we will use a value of 
a = 1/3. 

Trujillo et al. (2001b) noted that as the expo- 
sure depth increases and the outer radius there- 
fore increases, this definition of light concentration 
loses its significance and tends to a value of 1 for 
all galaxies. For the Sersic models, this concen- 
tration index can be written as a function of the 
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2n, bn [aj^^ 



Ijr. 



2n 



\ ' e,mod / 



1/n 



(5) 



2.3. A new concentration index 

We describe here a 'third galaxy concentration' 
(TGC) index defined in Trujillo et al. (2001b) as 



TGC{a) 



(6) 



Here, E{rg) means the isophote which encloses 
half of the total light of the galaxy^, and E{are) 
is the isophote at a radius a (0 < a < 1) times 
Vg. Again, this is a flux ratio. For a Sersic profile 
which extends to infinity, 



TGC{a) 



7(2n,^ai/") 
7(2n, bn) 



(7) 



For a range of different values of a, this concentra- 
tion index is shown in Figure 1 to increase mono- 
tonically with the value of n, revealing a relation 
between the central galaxy concentration of light 
and the global galaxy structure as defined by the 
shape parameter n (see Trujillo et al. 2001b). That 
is, the shape parameter n can be thought of as 
more than just a parameter that describes the cur- 
vature of the light profile, but has the additional 
physical meaning that it describes the degree of 
central concentration of at least luminous matter 
in Elliptical galaxies and bulges. This relation be- 
tween the central luminous concentration and the 
global structure of the galaxy can be made even 
more succinctly by deprojecting the various r^/" 
profile models to obtain their spatial luminosity 
density profiles (Ciotti 1991) for different values 
of n (Figure 2). One can immediately see the in- 
creasingly dramatic rise in central density^ with 
increasing n. 

For larger values of a (e.g. >0.5) the TGC index 
loses its ability to clearly distinguish galaxies with 
different structural profile shapes (i.e. different n). 



^In practice, Te is the observed half-light radius re,obs 
not r-e,mod- 

^For n > 1 the central density is actually infinite, occurring 
at a singularity (Ciotti 1991). 
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Low values for a, such as 0.2, provide a good range 
to the TGC index for different n, but in reahty 
are not so practical (especially when dealing with 
high-redshift galaxies). We therefore propose to 
use a = 1/3 (keeping the 3:1 ratio of the previous 
concentration indices). 

Equation 3 shows how reducing the outer ob- 
served radius r^^ of an r^/" profile reduces the 
observed galaxy luminosity. This in turn reduces 
the observed effective half-light radius from r,,.mod 
to re,obs (see Trujillo et al. 2001b, their Figure 9). 
As a result, the TGC index depends on the outer 
radius in the following way 
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[2"'^" (;:;:") 


1/n" 



where re,obs/''e,mod must first be computed given 
''fin/''e,mod (Trujillo et al. 2001b, their section 5.2). 

3. Stability Analysis of the VEirious con- 
centration indices 

In Figure 3 we have shown how the above three 
concentration indices vary as the outer final ra- 
dius is varied. One can see that the commonly 
used concentrations indices do not simply reveal 
the luminous structural concentration one hopes 
to measure but can be heavily biased by the ra- 
dial extent used to compute them. Consequently, 
concentration indices derived for exactly the same 
galaxies observed with first a 'shallow' and then 
a 'deep' exposure will be different. Clearly the 
"mean concentration index" (Section 2.2) is unre- 
liable by itself to provide any kind of meaningful 
galaxy classification. Okamura et al. (1984) rec- 
ognized this short-coming and, despite their hope 
to define a single fundamental structural quantity, 
had to resort to the introduction of an additional 
parameter, namely, the mean surface brightness 
of each galaxy. Even then they still had to ac- 
knowledge that the range of different global pro- 
file shapes (due also to the different bulge/disk 
combinations which galaxies possess) introduced 
additional scatter which left the previous two pa- 
rameters still unable to provide an accurate clas- 
sification of galaxy types along the morphologi- 
cal sequence (their Figure 3). They were however 
able to accurately (85%) classify galaxies as either 



early- type (E-SO/a) or late- type (Sb-Im). This 
broad categorical restriction was confirmed by Doi 
et al. (1992), Doi, Fukugita, & Okamura (1993), 
and then again by Abraham et al. (1994) who 
explicitly dealt with ellipticity and extended the 
formalism to the study of small and faint images 
where individual pixel information is important. 
Doi et al. (1993) also showed that the effects of 
seeing can result in one not even being able to dis- 
tinguish between these two broad categories when 
a small number of scale-lengths are sampled. Most 
recently, Bershady, Jangren, & Conselice (2000)^ 
have introduced a third parameter, namely color 
- which Hubble (1936) noted was correlated with 
morphological type - that has enabled them to 
separate galaxies into three classes: early, inter- 
mediate, and late-type. 

While the "mean concentration index" in com- 
bination with other galaxy parameters may be 
able to broadly categorize galaxies, it should 
clearly not be used as a tracer of the concentration 
of the luminous matter in galaxies. Its value de- 
pends sensitively on the depth of the image, and 
is therefore not closely related to any underlying 
physical property of the galaxy. Even for the same 
galaxy cluster, where the exposure details are the 
same, the situation is a mess. Figure 3 reveals 
that this index cannot even distinguish between 
a pure exponential disk (for e.g. an Sd galaxy) 
observed to three effective radii and a giant r^^* 
elliptical galaxy observed to one effective radii. 

The original concentration index defined as the 
ratio between radii enclosing different fractions (or 
quartiles) of the total galaxy luminosity (Fraser 
1972; dc Vaucouleurs & Aguero 1973; Fraser 1977; 
de Vaucouleurs 1977) is significantly more stable 
than the "mean concentration index" . This is 
because this index is less affected by where the 
galaxy profiles are taken to terminate. Although, 
it should be noted that this index is noticeably 
less well behaved when one removes the logarithm 
used in Figure 3. On the other hand, the TGC 
index appears to be very stable. 

In order to compare the C31 index with the 
TGC index, we have plotted in Figure 4 the radii 
at which the two indices get within 10% of their 

^Following Kent (1985), Bershady et al. (2000) used a ratio 
of radii for their concentration index C, such that C = 
51og[r(80%)/r(20%)]. 
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asymptotic value, which occurs at an infinite ra- 
dial extent. One can see that the TGC index ac- 
quires this level at a reasonable number of effective 
radii. The C31 index performs significantly worse. 
We therefore conclude that of the various concen- 
tration indices, the "mean concentration index" 
should not be used for the study of one-component 
systems (i.e. Elliptical galaxies, dwarf Elliptical 
galaxies, Spiral galaxy bulges, pure exponential 
disks). It's application to two-component systems 
may be even less appropriate and will be studied 
in a forth-coming paper. The TGC index should 
be used in preference to the C31 index. The rel- 
ative independence of the TGC index on the ob- 
served radial extent, and therefore on the individ- 
ual intrinsic galaxy surface brightness and expo- 
sure details, may mean that this single parameter 
can be used on its own to quantify, at least one- 
component systems, without the need to obtain 
and calibrate surface brightness levels and colors; 
that is, even non-photometric data could be ana- 
lyzed. 

3.1. Tests with observed galaxy profiles 

In the above section we have analyzed the sta- 
bility of the various concentration indices against 
truncations in the radial extent of model Sersic 
light-profiles. We now perform one of what is no 
doubt many possible tests which has no depen- 
dency on the Sersic model and uses real galaxy 
profile data containing noise, sky-subtraction er- 
rors, and possible deviations from perfect r^/" 
models. 

In practice, to determine the total galaxy light 
attempts arc made to account for the fraction of 
light which may reside outside the last measured 
isophote. In order to do this, the growth curve is 
normally extrapolated by one of a variety of tech- 
niques, such as a convenient mathematical func- 
tion, or fitting some ad-hoc model, or even some- 
times by eye. The extrapolation term, and the 
corresponding uncertainty, depend on a variety of 
factors, such as the shape of the light profile, the 
depth of the galaxian image, the accuracy of the 
sky-background subtraction, the particular way 
the extrapolation is computed, and so on. For this 
reason, it is of interest to compare (in a model- 
independent way) the measured concentration in- 
dices using observed profiles, and to check which 
index is least sensitive (i.e. more robust) to the 



above potential sources of error in estimating the 
total galaxy light. The "mean concentration in- 
dex" will not be computed here as it has already 
been shown to be a poor estimator of concentra- 
tion, equal to 1 for any profile measured to a large 
radial extent. 

We have analyzed the Virgo galaxy profiles pre- 
sented in Caon et al. (1993, 1994) through the fol- 
lowing experiment: 

1. The galaxy growth curve was computed, out 
to the outermost measured radius, using the 
observed major-axis light profile and the el- 
lipticity and position angle profiles. The to- 
tal B-band magnitude was derived by ex- 
trapolating the growth curve to infinity (see 
the above two references) to give jtib. 

2. Next, a variable magnitude Am was added 
to iriB to give mtot- Negative values of Am 
correspond to changes in the magnitude due 
to an overestimate of the total galaxy lumi- 
nosity; positive values may represent trunca- 
tions in the light profile, or represent magni- 
tudes within some isophotal threshold (such 
as V26), or simply account for errors from 
the true total galaxy magnitude. 

3. From the value of Wtot we then derived the 
effective half-light radius Tes and computed 
the TGC index. We also derived r25 and rj^ , 
the radii at which the enclosed luminosity is 
25% and 75% of the total luminosity, and 
determined the value of C31. Both concen- 
tration indices were derived independently 
of any profile model. 

4. Steps 2 and 3 were repeated for Am ranging 
from -0.2 to 0.4. 

By plotting these model-independently derived 
TGC and C31 indices as a function of Am, we 
were able to sec how they reacted to uncertain- 
ties/changes/errors in the total galaxy magnitude. 

For a given variation in galaxy profile shape, 
and therefore galaxy concentration, C31 may 
change its value by x% while the TGC index 
changes its value by y%, or vice- versa. There- 
fore the stability of different indices can not be 
compared in terms of their percentage changes, as 
different percentage changes may accurately re- 
flect identical changes in galaxy structure. The 
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sensitivity of each index to nicasurement errors 
should therefore perhaps be viewed in the Hght of 
the impUed changes in galaxy structure. Figure 5 
shows the derived TGC and C31 indices plotted 
against Am for a random sub-sample of galaxies 
possessing a range of Sersic indices n. It clearly 
shows that the TGC index is more stable than 
C31 because, for the same change in Am, the 
TGC index spans a smaller interval in An, that 
is, a smaller change in galaxy structure. In other 
words, if one used the model-independent deriva- 
tions for the TGC and C31 indices as a way to 
estimate the true galaxy structure as represented 
by n, the TGC index would give a more stable 
estimate which is less prone to uncertainties in 
rritot than the index C31 is. 

The above point is dramatically illustrated in 
Figure 6 and Figure 7. Figure 6 shows the TGC 
index derived from the best fitting Sersic model 
to the observed galaxy light profile and the TGC 
index derived directly from the light profile data 
itself, with no dependence on a r^/" model. In 
the model-independent case, the estimated to- 
tal galaxy magnitude was then increased and de- 
creased by Am=0.2 mag (i.e. spanning a range of 
0.4 mag) and the TGC index re-computed. The 
first thing one notes from Figure 6 is that the 
model-dependent and model-independent values 
agree reasonably well with each other. The second 
point, which receives emphasis when one simulta- 
neously considers Figure 7, is that the range in 
values for the TGC index is quite well constrained 
when one varies the total galaxy flux by nearly 
some 50%. The C31 index does not behave any- 
where near as well (Figure 7). It can be clearly 
seen to be far less stable to errors in the total 
galaxy flux than what the TGC index is. 

4. Correlations between galaxy concentra- 
tion, structure, and dynamics 

The concentration index, in its various guises, 
has been shown to correlate, albeit sometimes 
poorly, with galaxy morphological type; indeed, 
Doi et al. (1993) suggested it be used, together 
with the observed mean surface brightness, as a 
means to morphologically classify different galaxy 
types (sec also Bcrshady et al. 2000 and references 
within). We show for the first time in the follow- 
ing subsections that the TGC index, at least for 



the family of Elliptical galaxies, is strongly corre- 
lated with all the fundamental galaxy parameters. 
The concentration appears to not only reflect the 
general morphological structure, but is intimately 
related with the total luminosity, size, brightness, 
and central velocity dispersion of a galaxy. 

One may ask, "But isn't the TGC index merely 
another way of expressing the exponent n?" Or, is 
this central concentration index (a quantity which 
can be measured independently of any model or 
value of n) a fundamental quantity which is inti- 
mately linked with the nature and evolution of El- 
liptical galaxies. What are the fundamental quan- 
tities which should be plotted against each other 
to gain insight into the nature of Elliptical galax- 
ies? At this point we don't know. We there- 
fore show, for a sample of Elliptical galaxies, cor- 
relations between concentration and: luminosity, 
surface-brightness, scale-size, and velocity disper- 
sion. What past correlations with n actually mean 
are somewhat vague. Only when you understand 
what the correlations mean are you able to say 
something about the physical behavior of galaxies. 
What we shall see below is that the larger, more 
luminous and massive galaxies are more centrally 
concentrated. 

4.1. Luminosity 

The previous analysis, in particular Figure 1, 
reveals that the central luminous concentration in 
Elliptical galaxies, and also Spiral galaxy bulges, 
must be related to their global luminous struc- 
ture. This is because the more luminous galaxies 
and bulges are known to possess larger values of n 
(Caon et al. 1993; Young & Currie 1994, 1995; An- 
dredakis et al. 1995; Jerjen & Binggeli 1997; Gra- 
ham 2001) and consequently they must also have 
higher central concentrations of light than the less 
luminous bulges. Figure 8 shows that this is in- 
deed the case. Here, the TGC index has been plot- 
ted against the bulge luminosity from a sample of 
Virgo dwarf Elliptical galaxies (Jerjen, Binggeli, & 
Freeman 2000) and a sample of Virgo and Fornax 
early-type galaxies (Caon et al. 1993; D'Onofrio 
et al. 1994). The Pearson correlation coefficient 
between the TGC index and luminosity is —0.82. 
Given our new understanding of the relationship 
between the Sersic index n and the central concen- 
tration of luminous mass, the correlation between 
total luminosity and concentration is perhaps not 
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surprising. It has, however, to the best of our 
knowledge, never been shown before for a sam- 
ple of Elliptical galaxies. It can also be derived 
completely independent of any light-profile model 
(and hence value of n), although here we have used 
the TGC index from the best-fitting Sersic model 
because of the similarity seen in Figure 6 and be- 
cause we do not have the images for the dwarf 
galaxy data set. 

To calibrate the luminosity density profiles in 
Figure 2 requires some measure of the central in- 
tensity or central surface brightness /xq which is, 
for an r^/"* model, 7.67 mag brighter than the sur- 
face brightness /ie at one effective radius. Now, 
the Kormendy (1977) relation for Elliptical galax- 
ies tells us that /Xe oc 3 log re, and therefore^ the 
magnitude M oc — 2 log re . Thus, assuming r^/* 
profiles, the more luminous galaxies should pos- 
sess larger effective radii and fainter central in- 
tensities, which equivalently implies lower central 
luminosity densities. 

It must be stressed that we are not referring 
to the behavior of the luminosity density of the 
core within the central ~arcsecond as revealed 
with HST resolution (Rest et al. 2001, and refer- 
ences within), but to the properties derived from 
the global galaxy profile. Recent findings show 
that many Elliptical galaxies possess supermas- 
sive blackholcs at their centers (Kormendy & Geb- 
hardt 2001 and references within). Gravitational 
slingshots of stars which come to close to the 
central massive blackhole, or coalescing massive 
blackholcs from the progenitors of a merger may, 
to varying degrees, evacuate the core of a bulge 
and ther(^by reduce the original inner light profile 
(Ebisuzaki, Makino, & Okamura 1991; Makino & 
Ebisuzaki 1996; Quinlan 1996; Quinlan & Hern- 
quist 1997; Fabcr et al. 1997; Milosavljevic & Mcr- 
ritt 2001). If mergers involve strong gaseous dissi- 
pation and central starbursts this may also modify 
the nuclear profile (Mihos & Hernquist 1994), as 
do adiabatic black hole growth models (van der 
Marel 1999a,b, 2001). In this paper we are, how- 
ever, not talking about the very inner density of 
HST-resolved nuclear cusp slopes which may have 
been modified by the central blackhole, but are re- 
ferring to the global galaxy structure as seen with 
ground based resolution. 

^For an r^/* law, m = ^le — 5 log re — 3.388. 



The estimated masses of the central blackholcs 
in bulges have been shown to positively corre- 
late with the total luminosity of the host galaxy 
(Kormendy 1993; Kormendy & Richstone 1995), 
which would imply, assuming r^^^ profiles, that 
their masses are greater for galaxies with lower 
central luminosity densities (again, we are not 
referring to the modified HST resolved cusps). 
Why it should be that galaxies which are glob- 
ally more disperse have greater central blackhole 
masses must be explained if one is to assume that 
all Elliptical galaxies follow the r^/^ law. How- 
ever, this problem is quickly dismissed when one 
realises that it is a somewhat artificial problem 
which was created from the simplification of galax- 
ies through the assumption of r^/^ profiles. Not 
only do a range of light profile shapes exist, as do 
varying degrees of concentration, but the central 
(observed with groimd-based telescopes) surface 
brightness of all but the brightest Elliptical galax- 
ies actually brighten with increasing galaxy mag- 
nitude (Figure 9)"'' in contradiction to the predic- 
tion of the Kormendy relation when coupled with 
the assumption of r^^^ profiles. That is not to 
say that the Kormendy relation is wrong. Indeed, 
the brighter galaxy sample members in Figure 9 
roughly follow the Kormendy relation, having a 
slope of ^3 in the /ie logrg plane (Figure 10). (Al- 
though, the smaller and fainter galaxies do not fol- 
low the Kormendy relation, but there is certainly 
no suggestion that fi^ brightens with r,..) Figure 9 
reveals that the mass of the central supermassive 
blackholes are therefore positively correlated with 
the central luminosity density of the host galaxy 
as derived from the global luminosity profile, with 
the caveat that the very inner light profile cusps 
have likely been re-shaped by the central massive 
blackhole. 

We wanted to be confident that the trend seen 
in Figure 9 could not be explained by the influence 
of atmospheric seeing on what might in fact be 
r^/* law profiles which have smaller half-light radii 
the fainter the galaxy magnitudes arc. To explore 
this we convolved a series of r^/** profiles having a 
range of half-light radii from 5 to 50 arcsec with a 
Gaussian PSF having a FWHM of 2 arcsec, com- 
parable to the worst seeing conditions under which 

^ There is some suggestion of a turn around at Bt < 10 in 
Figure 9, evident amongst the brightest galaxy members in 
Figure 6c of Faber et al. 1997. 
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the central CCD data was obtained for the galaxy 
sample shown in Figure 9. We found that the av- 
erage surface brightness within the inner circle of 
radius 1 arcsec was under-estimated by ^0.5 mag 
arcsec"^ when re=5" and by ~0.3 mag arcsec"^ 
when re=50". These differences where of course 
even smaller when the values within the inner 2" 
radius were considered, and the trend seen in Fig- 
ure 9 did not vary noticeably when we used the 
observed data within the inner 2". The effects of 
seeing are therefore unable to resurrect the possi- 
bility that the trend seen in Figure 9 is compatible 
with all galaxies having r^/* law profiles. 

4.2. Central surface brightness and effec- 
tive radii 

Khosroshahi, Wadadekar, & Kembhavi (2000) 
and Mollenhoff & Heidt (2000) have modelled 
the light profiles from a sample of spiral galaxies 
with a seeing-convolved r^^"^ bulge and exponen- 
tial disk model. The strongest correlation found 
between any of the structural parameters explored 
by these authors was between the central bulge 
surface brightness and the bulge shape parameter 
n. Khosroshahi et al. obtained a linear correlation 
coefficient r = —0.88 at a significance level bet- 
ter than 99.99%, and Mollenhoff et al. obtained a 
value of r = —0.86. Positive correlations between 
n and were also found by these authors. Thus, 
from Section 2, the central luminosity density of 
spiral galaxy bulges (as represented by /xq) must 
be positively correlated with not only n but also 
with the bulge concentration of luminous matter. 
One may indeed ask which is the more fundamen- 
tal connection. 

We show in Figure 11 and Figure 12 the cor- 
relation between the TGC index and re and /zq 
for our sample of dwarf Elliptical and Elliptical 
galaxies. The behavior in Figure 11 is similar to 
that already known between n and rg (e.g. Caon 
et al. 1993; Graham et al. 1996). The correlation 
coefficient is 0.67, and turns out to be the weakest 
correlation we find between concentration and any 
of the other galaxy parameters presented here. 

In Trujillo et al. (2001b) we showed that the 
range of different galaxy structures which exists 
amongst the Elliptical population cannot be due 
to parameter coupling in the r^^" models, and sys- 
tematically varies with model-independent quan- 
tities such as effective half-light radius and lumi- 



nosity. We caution, however, that for large values 
of n the effects of seeing on the light profile can be 
substantial at small radii (Trujillo et al. 2001a). 
We also note that while the r^/" profiles arc good 
at describing the global profile shape, they can re- 
quire modification at some small inner radii (Jaffe 
et al. 1994; Fcrrarcscct al. 1994; Laucr et al. 1995; 
Faber et al. 1997). The notably bright (~< 15 
mag arcsec"^) central (r=0) surface brightnesses 
expected from the very high values of n are per- 
haps also unlikely to be realised. Either the Sersic 
model is no longer appropriate to describe the very 
inner profile, and/or the extremely high inner den- 
sities modify the actual profile. The presence of 
(and likely past infall of material into) a central 
supermassivc BH is also likely to disturb the cen- 
tral cusp slope, as can mergers and other mech- 
anisms (Makino & Ebisuzaki 1996; Faber et al. 
1997; Quinlan & Hcrnquist 1997; Merritt & Quin- 
lan 1998). We therefore note that Figure 12 (see 
also Jerjen & Binggeli 1997), showing the relation 
between the TGC index and the central galaxy 
surface brightness obtained from the best-fitting 
r^/" model, should be regarded as somewhat pre- 
liminary and a more detailed investigation shall 
be forthcoming. Despite these words of caution, 
excluding those galaxies with values of n greater 
than 4, that is, removing those galaxies which may 
have overly bright central surface brightness esti- 
mates, still resulted in a correlation coefficient of 
-0.88 between the central concentration and cen- 
tral surface brightness. (Using all of the galax- 
ies gave a correlation coefficient r = —0.94.) We 
stress again that if all Elliptical galaxies followed 
the r^l'^ law, then the concentration index would 
be the same for all galaxies, and no correlation 
woiild exist between concentration and any of the 
other galaxy parameters; this is similarly the case 
if all dwarf Elliptical galaxies and Spiral galaxy 
bulges were to possess the same universal profile. 

4.3. Velocity dispersion and mass 

Combining the luminosity - central concentra- 
tion relation (Figure 8) with the Faber-Jackson 
(1976) relation between luminosity and central ve- 
locity dispersion immediately implies that the cen- 
tral velocity dispersion and therefore mass - as far 
as the central velocity dispersion is a measure of 
the galaxy mass - must positively correlate with 
the galaxy concentration. Using the early-type 
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galaxies from Caon et al. (1993) and D'Onofrio ct 
al. (1994) which have available central velocity dis- 
persion measurements from Hyper cat®, Figure 13 
shows that the TGC index is strongly (r = 0.80 
for the Ellipticals) correlated with the velocity dis- 
persion, even with the heterogeneous nature of the 
dynamical data. Excluding galaxies with central 
velocity dispersions less than 100 km s~^ did not 
change things appreciably, nor did using the cen- 
tral velocity dispersion catalog of McElroy (1995) 
where we found a value of r = 0.82. This is a fun- 
damental result; together with the previous corre- 
lations it tells us that the more massive a galaxy is, 
the more centrally concentrated it must be. The 
virial theorem (and its observational counterpart 
the Fundamental Plane: Djorgovski & Davis 1987; 
Dressier et al. 1987) relates the luminosity, size, 
and velocity dispersion (kinetic energy) terms, but 
does not explain why greater mass should imply 
higher central concentration. Theories of gravita- 
tional collapse and galaxy formation must be able 
to explain this. 

We speculate here that galaxy concentration 
may provide the physical mechanism for the ob- 
served connection between the stellar velocity dis- 
persion (mass) of a bulge and the mass of its cen- 
tral supermassive BH (Merritt 2000; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000; Merritt & Fer- 
rarese 2001; Marconi et al. 2001; Sarzi et al. 2001). 
The exact process which explains the reason for 
the existence of this correlation is not yet known, 
although many theories have been proposed (Ef- 
stathiou & Rees 1988; Haehnelt & Rees 1993; 
Ciotti & Ostriker 1997; 2001; Haiman & Loeb 
1998; Silk & Rees 1998; Blandford 1999; Haehnelt 
& Kauffman 2000; Kauffman & Haehnelt 2000; 
Ostriker 2000; Adams, Graff, & Richstone 2001). 
Bigger galaxies with higher central concentration 
and stronger potential wells would naturally sup- 
ply more fuel to their inner regions. This could 
take the form of more efficiently funneling gas to 
build the central accretion disks that likely feed 
the quasars wc; obsc; it high-redshift (Soltan 
1982), perhaps leaving remnants such as the nu- 
clear disks we observe today (Rest et al. 2001 and 
references within) which would now encircle the 
inactive heart of what was once a quasar. 

'Hypercat can be reached at http://www-obs.univ- 
lyonl.fr/hypercat/. 



It is not clear why simply having a higher veloc- 
ity dispersion alone can be the fundamental physi- 
cal mechanism for greater BH masses. On the one 
hand it does imply that more material (mass) is 
in the galaxy to build/feed a blackhole, but we 
saw in Section 4.1 that if all galaxies were de- 
scribed by the r^^^ law then galaxies with greater 
mass (higher velocity dispersion) would have lower 
central luminosity densities and be more disperse. 
This would require a low-density environment to 
either favor the formation of more massive black- 
holes or to be the product of their evolution. 
(Once again, we stress that we arc not referring to 
the very inner nucleus as revealed with HST imag- 
ing.) We propose here that the presence of mass 
itself may not be the end of the story. Figure 13 
shows that more massive galaxies are more cen- 
trally concentrated than less massive galaxies. We 
suggest that exactly how the galaxy mass is dis- 
tributed/concentrated may be an important fac- 
tor. This could be tested through a photomet- 
ric campaign which measures the concentration in 
those galaxies with blackhole mass estimates, and 
determines the strength and scatter of the corre- 
lation between these two quantities. Although, 
Ferrarese & Merritt (2001) and Gebhardt et al. 
(2000) found that the intrinsic scatter in the BH 
mass velocity dispersion diagram is small or neg- 
ligible (i.e. consistent with the measurement er- 
rors alone) for those galaxies with most reliable 
blackhole mass measurements. This would then 
leave no room for improvement with a BH mass - 
galaxy concentration diagram, and would suggest 
that the former is indeed the more fundamental re- 
lation. This result may hold firm, although it will 
be of interest to see what happens when more BH 
mass estimates are obtained and refined on sev- 
eral fronts. Improvements will come with the use 
of non-axisymmetric dynamical models, improved 
knowledge of line broadening mechanisms (Barth 
et al. 2001), and addressing concerns that the dy- 
namical effects of some blackholes may occur at 
resolutions lower than presently probed (Qian et 
al. 1995; de Zeeuw 2000). Additional corrections 
for the finite slit width of STIS on the HST (Ma- 
ciejewski & Binney 2001) may also prove crucial, 
and significantly lower the current blackhole mass 
estimates. 
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5. Conclusions 

We have explored the stabiUty of several differ- 
ent galaxy light concentration indices as a function 
of galaxy exposure depth, or rather, the number 
of effective radii sampled. This analysis has been 
confined to one-component stellar systems such as 
normal Elliptical galaxies, dwarf Elliptical galax- 
ies, Spiral galaxy bulges, and exponential disks. 
Our investigation reveals that the "mean concen- 
tration index" , often used for parametrizing faint 
and high-redshift galaxies, is a horrendously poor 
estimator of galaxy light concentration and it's use 
for such a task on its own should be abandoned. 
To illustrate this claim, we have shown that this 
index is unable to distinguish between a giant r^l^ 
Elliptical galaxy measured to one effective radii 
and a pure exponential disk measured to three ef- 
fective radii. The de Vaucouleurs C31 index per- 
forms notably better but is still heavily dependent 
on the outer galaxy radius one reaches. The cen- 
tral concentration index introduced in Trujillo et 
al. (2001b) is shown to be the more stable of the 
indices, changing in value (with increasing galaxy 
radius) by less than 10% once a few effective radii 
have been sampled, and is more robust against 
measurement errors. 

Given that Elliptical galaxies and bulges are 
not homologous systems, they therefore possess 
a range of different (light) concentrations. The 
global profile shape, which can be parameterized 
by the value n from the Sersic r^/" law, is inti- 
mately connected with the degree of galaxy light 
(and mass) concentration. Which of these two 
quantities is the more fundamental is not clear; the 
latter quantity can however be measured directly 
from the image or light-profile, derived indepen- 
dently of any galaxy model. For a sample of dwarf 
Elliptical galaxies and normal Elliptical galaxies 
we have presented strong correlations between the 
central concentration index, as defined in this pa- 
per, and the global galaxy parameters: luminos- 
ity, effective radius, and central surface brightness. 
We also present the first ever correlation between 
galaxy concentration and velocity dispersion for 
a sample of Elliptical galaxies, showing that the 
more massive galaxies arc more centrally concen- 
trated; this should provide a valuable clue into the 
physics of gravitational collapse and galaxy forma- 
tion. Lastly, we speculate and provide a means to 



test that the central concentration of at least lu- 
minous matter in Elliptical galaxies, that is, how 
this matter is distributed, may be an important 
quantity regarding the formation of supermassive 
BHs. 

We arc happy to thank Helmut Jcrjen for pro- 
viding us with the dwarf Elliptical data used in 
Figure 8-12. We also wish to acknowledge and 
thank the anonymous referee for their suggestions 
and comments which helped to improve this pa- 
per. 
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Fig. 1. — The central galaxy concentration index 
as defined in equation 6 and equation 7 is shown 
as a function of the Sersic shape parameter n, 
for different values of a. (Extension of Figure 3 
from Trujillo et al. (2001b) who plotted values of 
a equal to 0.3 and 0.5.) 



Fig. 2. The dcprojcctcd (spatial) luminosity- 
density profiles for a range of Sersic r^^" mod- 
els are shown as a function of the deprojected 
radius s, normalized at 1 from the projected 
r^/" model. The increased central concentration 
of these models with n is strong and obvious. 
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Fig. 3. — Various concentration indices have been 
plotted, for different galaxy profile shapes n, as a 
function of the outer radius rfjn used to compute 
each index (see text for details). Panel a) presents 
the concentration index as defined by Abraham 
et al. (1994). Panel b) presents the logarithm of 
the concentration index defined be de Vaucouleurs 
(1977). Panel c) presents the concentration index 
defined by Trujillo et al. (2001b). Panel d) to f) 
show the various indices as a function of the ob- 
served half-light radius, which changes with rfin- 
The value of a used to compute C31 and TGC has 
been set to 1/3. The TGC index is clearly the 
most stable of the three, while the index in panel 
a) is horrendously unstable. 
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Fig. 4. — Stability analysis of the the C31 concen- 
tration index and the TGC index. The radius rgn 
(normalized by the (effective radius rc.mod of each 
^i/n jjjodel) where the concentration indices get 
within 10% of their maximum value (which occurs 
at an infinite radial extent) is shown as a function 
of n. The TGC index is clearly more stable with 
radius, and hence exposure depth, than the C31 
index. 
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Fig. 5. — TGC and C31 indices (computed inde- 
pendently of any Sersic r^/" model) for a sample of 
real galaxy profiles spanning a range of structural 
shapes (i.e. n). The light-profile data was taken 
from Caon et al. (1990) and (1994) for NGC 4623, 
NGC 1379, NGC 4636, and NGC 4365, in order of 
increasing galaxy light concentration. The influ- 
ence on the concentration indices by supplement- 
ing the observed galaxy magnitude within some 
outer radius with Am is revealed (see text for de- 
tails). The term Am can represent the magnitude 
beyond the outer radius which was missed, or rep- 
resent truncations to an inner isophotal level, or 
one of a number of sources of observational er- 
ror. Dotted horizontal lines mark the values of 
the TGC and C31 indices corresponding to infi- 
nite Sersic models with n = 0.5, 1, 2. . .10. 
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Fig. 6.— The stabihty of the TGC index has 
been explored with real galaxy data by deriving 
its (model-independent) value when an error of 
±0.2 mag is added to the (model-independent) 
estimated total galaxy magnitude. The x-axis is 
the TGC index derived from the best-fitting Sersic 
model. Only galaxies with major-axis fits which 
were marked as 'Good' (that is, not those whose 
quality was marked as 'Fair' or 'Poor') from Table 
2 of Caon et al. (1993) have been used. 
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Fig. 7. — Same as Figure 6 but for the C31 index. 
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Fig. 8. The TGC index is plotted against galaxy 
magnitude for a sample of Virgo dwarf Elliptical 
galaxies (triangles) taken from Jerjen, Binggeli, & 
Freeman (2000) and a sample of early-type Virgo 
and Fornax galaxies (circles) taken from Caon 
et al. (1993) and D'Onofrio et al. (1994). The 
TGC index has been derived using equation 7, and 
the magnitudes are the observed (i.e. not model- 
dependent) values. All galaxies which could be 
modelled with a Sersic profile along the major-axis 
have been included. 



Fig. 9. — The observed (i.e. model-independent), 
central mean surface brightness within a circle of 
radius 1 arcsec is plotted against the observed 
(model-independent) galaxy magnitude. Filled 
circles represent Elliptical galaxies and open cir- 
cles represent SO galaxies from the complete Virgo 
and Fornax galaxy sample of Caon et al. (1993) 
and D'Onofrio et al. (1994) observed from the 
ground. 
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Fig. 10. — The (model-independent) surface 
brightness fie at the (model-independent) effec- 
tive half-light radius re is plotted against fg. The 
data are from the complete Virgo and Fornax 
galaxy sample studied in Caon et al. (1993) and 
D'Onofrio et al. (1994). 
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Fig. 12. — The TGC index shown in Figure 8 has 
been plotted against the central surface brightness 
of each galaxy. The only two galaxies with val- 
ues for n greater than 10 have been excluded (see 
text for additional comments). Circles represent 
the early-type galaxies, and triangles represent the 
dwarf Elliptical galaxy sample. 



Fig. 11. — The TGC index shown in Figure 8 is 
plotted against each galaxy's (model-independent) 
effective radii. Circles represent the early-type 
galaxies, and triangles represent the dwarf Ellipti- 
cal galaxy sample. 
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Fig. 13.— The TGC index shown in Figure 8 
has been plotted against the central velocity dis- 
persion for the early-type galaxies with published 
kincuiatical data. Filled circles represent the El- 
liptical galaxies, and open circles represent SO 
galaxies. 
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